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Abstract 

A recent study on axisymmetric ideal magnetohydrodynamic equilibria with 
incompressible flows [H. Tasso and G. N. Throumoulopoulos, Phys. Plasmas 
5, 2378 (1998)] is extended to the generic case of helically symmetric equi- 
libria with incompressible flows. It is shown that the equilibrium states of 
the system under consideration are governed by an elliptic partial differential 
equation for the helical magnetic flux function ip containing five surface quan- 
tities along with a relation for the pressure. The above mentioned equation 
can be transformed to one possessing differential part identical in form to 
the corresponding static equilibrium equation, which is amenable to several 
classes of analytic solutions. In particular, equilibria with electric fields per- 
pendicular to the magnetic surfaces and non-constant-Mach-number flows are 
constructed. Unlike the case in axisymmetric equilibria with isothermal mag- 
netic surfaces, helically symmetric T = T{ip) equilibria are over-determined, 
i.e., in this case the equilibrium equations reduce to a set of eight ordinary 
differential equations with seven surface quantities. In addition, it is proved 
the non-existence of incompressible helically symmetric equilibria with (a) 
purely helical flows (b) non-parallel flows with isothermal magnetic surfaces 
and the magnetic field modulus being a surface quantity (omnigenous equi- 
libria). 
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I. Introduction 



Understanding the equilibrium properties of magnetically confined plasmas 
is one of the basic objectives of fusion research (see citations of Ref. JI|). In 
particular, in the last fifteen years, equilibria with flows which are induced 
when either neutral beams or/and electromagnetic power are employed to 
heat the plasma of magnetic confinement systems, has become an issue of 
increasing interest. These flows are usually associated with radial electric 
fields which play a role in the transitions to the improved confinement regimes 

I, 111. 

In a series of papers [7], |3| we investigated the magnetohydrody- 

namic (MHD) equilibrium of ideal plasmas with incompressible flows and 
translational as well as axial symmetry. The main conclusions of the above 
mentioned studies are the following: (a) If the equilibrium flows of cylindri- 
cal plasmas with arbitrary cross- sectional shapes are purely poloidal, they 
should be incompressible, (b) Exact equilibria were constructed for con- 
stant poloidal flow Mach numbers M 2 (see definition (|24l) in Sec. II) for 
both cylindrical and axisymmetric configurations, (c) For the physically ap- 
pealing class of cylindrical equilibria with isothermal magnetic surfaces their 
cross-sections must be circular at large, while no restriction is imposed on 
the magnetic surfaces of axisymmetric T = T(ip) equilibria, apart from the 
vicinity of the magnetic axis where they should be circular. 

The aim of the present report is to extend our previous studies to the 
generic case of helically symmetric ideal MHD equilibria with incompressible 
flows which describe the steady states of the plasmas of straight stellarators. 
It is shown that the equilibrium states satisfy an elliptic partial differen- 
tial equation for the helical magnetic flux function ip containing five surface 
quantities along with a relation for the pressure. This is the object of Sec. 

II. In Sec. Ill the above mentioned differential equation is transformed to 
one possessing differential part identical in form to the corresponding static 
equilibrium equation, which permits the derivation of several classes of ana- 
lytic solutions with non- constant M 2 and differentially varying electric fields 
perpendicular to the magnetic surfaces. We obtain exact solutions with flows 
in the symmetry direction, parallel to the magnetic field and in arbitrary di- 
rection. Equilibria with isothermal magnetic surfaces are then examined in 
Sec. IV. Sec. V summarizes our conclusions. 
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II. Equilibrium equations 



The ideal MHD equilibrium states of plasma flows are governed by the 
following set of equations, written in standard notations and convenient units: 

V ■ (pv) = (1) 

p(v • V)v = j x B — VP (2) 

V x E = (3) 
VxB=j (4) 
V-B = (5) 

E + vxB = 0. (6) 

The system under consideration is a helically symmetric magnetically con- 
fined plasma with flow. To describe this configuration we introduce cylin- 
drical coordinates R, z, with z along the rotation axis; helical symmetry 
implies that any physical quantity depends only on R and u = l<ft + kz. 
Though a helix is characterized by just one parameter, for instance a = k/l, 
it is convenient to keep both / and k in order to discuss the two limiting cases 
/ = corresponding to axisymmetry and k = corresponding to translational 
symmetry. Also, the vector 

= Nz-kR*V<p 

P + k 2 R 2 [) 

is introduced which is tangent to the helix r = const, u = const., with 
Vm = IVcf) + kVz. It then follows from Eq. (0) that h ■ S7G(R, u) = for 
any function G(R,u). Also, h has the properties 

—Oh] 

V'h = 0. Vxh^^.^h. (8) 

The divergence-free fields, i.e. the magnetic field B, the current density j and 
the mass flow pv can be expressed in terms of the stream functions ip(R, u), 
I(R,u), F(R,u) and Q(R,u) as 

B = /h+hxVV>, (9) 
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j = {£4 - 2klh 2 I)h - h x V/ 



and 



pv = 0h + h x VF, 
where £ is an elliptic operator defined by 



(10) 
(11) 
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;i2) 



and h 2 = h • h = I /{I 2 + k 2 R 2 ). Eqs. ([j])-® can be reduced by means 
of certain integrals of the system, which are shown to be surface quantities. 
To identify three of these quantities, the time independent electric field is 
expressed by E = — V$ and the Ohm's law @ is projected along h, B and 
VV>! respectively, yielding F = F(ip), $ = $(■?/>) and 
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[IF' - 0) = 



(13) 



(The prime denotes differentiation with respect to tp). A fourth surface 
quantity is derived from the component of momentum conservation equation 
(g) along h: 

(F') 2 \ 



I I 
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h 2 



(14) 



From Eq. (14) it follows that, unlike the case in static equilibria, I is not 
a surface quantity. This implies that in flowing plasmas the current density 
does not lie on magnetic surfaces [see Eq. ([TO])]. Solving Eqs. ( p~3|) and ([TJ]) 
for / and one obtains 

_ X - F'^'jh 2 



and 



Also, on account of Eq. 
form 







1 - {F'fjp 

F'X - p$'/h 2 
1 " {F') 2 /P ' 



(15) 



(16) 



I3|) , the velocity [Eq. ([□])] can be written in the 



F' 



B 



h 2 



(17) 
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With the aid of Eqs. (|TT|)-(]T7|), the components of Eq. (|2]) along B and 
perpendicular to a magnetic surface are put in the respective forms 



B 



v 
~2 







VP 







(18) 



and 



1 - 



+ 



r\2- 



+ 



F'F" 



P 



h 2 \Vij\ 2 -2klh 4 X\ |V# 



t (vv 2 - h 2 v(e/ P ) 2 ) + h 2 ^f + vp 



vv> = o 



(19) 



It is pointed out that Eqs. (p~8"D and (p~9[ ) are valid for any equation of state 
for the plasma. 

In order to reduce the equilibrium equations further, we employ the in- 
compressibility condition 

V ■ v = 0, (20) 
which on account of V • pv = implies that the density is a surface quantity 



P = /#)• 



(21) 



Consequently from Eqs. (|T3j) and flT4]) it follows that, similar to axisymmet- 
ric incompressible equilibria with purely poloidal flows [|l|] , helically symmet- 
ric plasmas with purely helical flows (0 = 0) can not exist; the only possible 
equilibria of this kind are cylindrical configurations with arbitrary cross sec- 
tional shapes 0. With the aid of Eq. (^1|), Eq. ([18]) can be integrated to 
yield an expression for the pressure, i.e. 



(22) 



We note here that, unlike in static equilibria, in the presence of flow magnetic 
surfaces in general do not coincide with isobaric surfaces because Eq. (Q) 
implies that B ■ VP in general differs from zero. In this respect, the term 
P s (i/>) is the static part of the pressure which does not vanish when v = 0. 
Eq. (|19D has a singularity when 



en 

p 



1. 



(23) 
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On the basis of the definitions v\ = — - fe_LY^j f or the helical velocity 

component, v\ h = — 1 ^ 1 for the Alfven velocity associated with the helical 
magnetic field, and the Mach number 



M 2 = ^- 

VAh 



(F 



r\2 



(24) 



Eq. (]23| ) can be written as M 2 = 1. If it is now assumed that ^ ^ 7^ 1 and 



Eq. (|22| ) is inserted into Eq. fll^D, the latter reduces to the elliptic differential 
equation 



1 - M 2 )C4> - -(M 2 )'|V^| 2 - 2klh 2 X 



+ 2 (l-M 2 ) + /^ ( Ps 



XF 1 ® 
1 — M 2 



1 / P(^) 2 ' 
2/i 4 VI — M 2 



0. 



(25) 



This is the equilibrium equation for a helically symmetric plasma with in- 
compressible flows. Eq. (p5|) contains the arbitrary surface quantities F(i[>), 
$(?/>), X(tJ)), p(%p) and P s {ip) which must be found from other physical con- 
siderations. Inspection of Eqs. (f25"D, (|T0| ) and ([Lj) shows that the singularity 
M 2 = 1 is the limit at which the confinement can be assured by the toroidal 
current density Ch alone. For M 2 7^ 1 the derivative of X 2 / (1 — M 2 ), related 
to the derivative of the magnetic field component in the symmetry direction 



by Eq. (13), partly compensates for the pressure gradient and inertial flow 
forces. 

It should be noted here that the equilibrium of helically symmetric equi- 
libria with incompressible flows was investigated in Refs. || [TO, ITTJ and [ITJ 



where a set of an elliptic differential equation and a nonlinear algebraic equa- 



tion was obtained (e. g. Eqs. (38) and (40) of Ref. 0]). Those equations, 
however, are coupled (and consequently one has to solve the one equation 
and then examine a posteriori whether the solutions are compatible with the 
other), a property which makes the derivation of analytic solutions tedious. 
On the other side, as concerns the present study, one has to solve a single 
differential equation [Eq. (p5|) 1 which is decoupled from Eq. ([22]). Once the 
solutions of Eq. ( ^5] ) are known, Eq. ( [Z2"D determines only the pressure. 
Moreover, as shown in section III, Eq. (|25"D can be transformed to one which 
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is amenable to several classes of analytic solutions describing equilibria with 
non-constant-poloidal-flow Mach numbers [(d/dU)M 2 ^ 0)]. 



III. Analytic equilibrium solutions 



Eq. fl25|) under the transformation [IJ|, [14J 

U{$) = [* [1 - M 2 (?//)] 1/2 M 2 < 1 



(26) 



reduces (after dividing by (1 — M 2 ) 1 / 2 ) to 



CU - Iklh 2 



X 



+ 



1 d 



X 2 



| ,11 P _ x dFd$ 
M 2 ) 1 / 2 ' 2dU \l-M 2 ) h 2 dU V s dtfdtf, 



1 d 
2/^df7 



/d$ N 

P U7, 



It is noted here that the requirement M 2 < 1 in transformation 

t2 tttVi rwr\ n t ( f\i / rC\ 1/2 



0. (27) 



ra) implies 



that t> 2 < v 2 , where v s = ( r yP/p) 1/2 is the sound speed. This follows from 



Eqs. 



and (in Gaussian units) 

(-y 



(7/2) 



8vrP 



1. 



Since, according to experimental evidence in stellarators ||, f|, |j, the (max- 
imum) value of the ion poloidal velocity (corresponding to Vh) in the edge 
region during the L-H transition is of the order of 1 Km/sec and the ion 
temperature is of the order of 100 eV, the scaling Vh <C v s is satisfied in this 
region. Therefore, the restriction M 2 < 1 is of non-operational relevance. 
Also, solving Eq. ([27|) only is sufficient to determine the equilibrium because 
it is not rj)(U) itself but Vip(U) = (1 - M 2 (U)) 1/2 VU which is needed to 
obtain any equilibrium quantity. 

Three classes of exact equilibria of Eq. (|2"T| ) can be constructed as follows. 

(a) Flows in the symmetry direction 
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This kind of equilibria correspond to F' — 0. From Eq. ( Jl5| ) it then follows 
that the function / becomes a surface quantity, / = X/(l — M 2 ), and Eq. 
fl2"?D reduces to 



£f/ - 2klh 




1 dP 



1 d 



M 2 j h 2 dU 

2" 



0. 



(28) 



Eq. (|28|) can be linearized and solved analytically for several choices of 
the surface quantities it contains. For example, two classes of polynomial 
solutions with respect to P and u can be constructed for (d/dU)[X(l — 



M 2 ) 



-1/21 



constant, (d/dU)P s = constant and (i) p(d&/dU) = const, and 



(ii) p(d$ / dU) 2 oc U. In case (ii), owing to the flow term containing the factor 
l/h 4 in Eq. (p7|) , the polynomials can describe configurations with either a 
single magnetic axis or two magnetic axes. A simple equilibrium with single 
magnetic axis corresponding to (d/dU)[X(l — M 2 ) -1 / 2 ] = is given by 



U 



R 2 (9R 



Q Z?2 p2 



P 4 - 5R\ 2 ) 



(29) 



where U c is the value of the flux function at the position of the magnetic 
axis (u = 0, P = P c ) and S is a parameter related to the shape of the flux 
surfaces. 



(b) Flows parallel to B 

Equilibria with B-aligned flows correspond to $' = 0. Eq. (^7|) then 
reduces to 



CU - 2Hh 2 



X 



1 d 



X 2 



(1 -M 2 )V2 2dU ll-M 2 



l_dP 
h 2 dU 



0. 



(30) 



which is similar in form to the equation governing static equilibria. (Analytic 
static helically symmetric equilibria have been derived in Refs. JL5|, [16|, |17|). 
Eq. (0) can be linearized and solved for (a) (d/dU)[X(l - M 2 )- 1 / 2 } = 
c + c\U and (d/dU)P s = d + diU, where c , c\, d and di are constant 
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quantities. The simplest solution corresponding to cq = c\ = d\ = is given 
by 

U = U c ^(2R 2 c -R 2 -5u 2 ). (31) 



For I = and k — 1 Eq. (plf ) reduces to the Hill's vortex axisymmetric 
equilibrium [pf|. 

(cj Non-parallel flows 

Equilibria of this kind are of particular interest because non-parallel flows 
with non-vanishing poloidal (helical) components are associated with radial 
electric fields which play a role in the transitions to improved confinement 
regimes in tokamaks and stellarators |4[]. Analytic, smooth equilibrium 
solutions of Eq. P7| ) can be derived for (i) p(d&/dll) 2 = constant and 
(ii) p(d§/dU) 2 = c 2 U, c 2 = constant. Because of the departure of the 
isobaric surfaces and current density surfaces from the magnetic surfaces, 
however, in case (i) it is not possible to obtain physically plausible equilibria, 
viz. equilibria with singly peaked density and pressure profiles (e. g. on 
the plane u = 0) vanishing on the plasma surface, a current density profile 
not possessing a component along V£7 on the plasma surface and finite E 
throughout the helical cross-section. In case (ii) a class of physically well 
behaved equilibria is derived with the ansatz 

X n dFd® TT 

Co = const, F s — X ———— oc U, 



(1 -M 2 )V2 ^ u ^ uu ' s dUdU 
M 2 = dU^ (ci = const., /I > 1) and p oc U x (0 < A < 1). 



Equation (|27|) then becomes 



CU - 2klh 2 c + + (32) 



A particular solution of Eq. (|3^) is given by Eq. (|29|) . The radial electric 



field on the plane u = following from Eq. (E^) is given by 



= _[/d-A)/2^ V£ (33) 

u=0 



Eq. (^) implies that |E#| vanishes at the magnetic axis R = Rc and on 
the plasma surface R = R s = [3(v^5 — 1) /2] l l 2 ]R c (which is determined by 



8 



U(u — 0, R) = 0) and has an extremum in between R c and R s , the position 
of which is a function of the parameters R c and A. This profile is consistent 
with those observed during the L-H transition in stellarators f|, |5J. 

III. Equilibria with isothermal magnetic surfaces 

For fusion plasmas the thermal conduction along B is fast compared to the 
heat transport perpendicular to a magnetic surface and therefore equilibria 
with isothermal magnetic surfaces are of particular interest. It is recalled that 
for cylindrical plasmas the relation T = T(ip) restricts the cross sections of 
the magnetic surfaces to be circular at large, while for axisymmetric plasmas 
no restriction is imposed on the shape of the magnetic surfaces apart from 
the vicinity of the magnetic axis where their cross sections should be circular. 

Under the assumption that the plasma obeys the ideal gas low P = RpT, 
Eqs. (0), (0) and ( [P7D lead to the following expression for the magnetic 
field modulus: 

|B| 2 = m + , (34) 

where E(ip) = 2(P S - P)p/(F') 2 and H(ip) = (p§'/F') 2 . Consequently, apart 
from the case of field aligned flows (H = 0), omnigenous equilibria, viz. 
equilibria with |B| being a surface quantity, are not possible. It is noted that 
this property also holds for axisymmetric equilibria ]T|]. 

Solving the set of equations (p5|) and ( |34] ) for |V?/>| 2 and Cip one obtains 



i\7/i* (wW i f^Y i(-a\u.iw -u k w\ 



(35) 



and 



C^-K^W-m-^-tf, (36) 

where i(if)), kijp), w(if)), fijp), g(i/i) and d(if)) are known functions of 

■0. With the introduction of the quantities x = R 2 , p = dip/dx, q = dip/du, 
r = d 2 ip/dx 2 and t = d 2 ip/du 2 , Eqs. (|35| ) and (^B) are written in the 
respective forms 

ixJ,2 + ^ q2 = 2{i+ h 3 + h k) (37) 

and 

Ax 1 111 

- r + il * p+ — t = - w --f-- g -- d . (38) 
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To integrate Eqs. (|37| ) and fl3"8| ) we apply a procedure suggested by Palumbo 
. Accordingly, employing i? and as independent coordinates instead of 
R and u (then w = u(x, ■0)), we have 



(39) 



and 



identically), Eq. ( p8[) reduces to 



dp 




<9p 




dp 


dx 






+ p 

</> 


dtp x 




dg 




dq 




t = 


du 




^ dip 


X 




3), 


(ID) and f + i' 



(40) 



Ax dp 
h 2 dx 



+ 4/ 2 p = -w - !(/ + - l(p + /) - l(d + m'). (41) 



Eq. (f4T|) can be considered as an ordinary linear first-order equation for 
p(x,ip) with respect to x and constant ip. Its solution is 



V 



1 . j2 . T(lp) I 2 + k 2 X 

'w+l 2 r(ip)) + 



Ak 2 x 



4 



4 



(g + j') + l 2 (d + m') + k 2 (d + m')^ 



where r(ip) is an integration "constant". 

Since u is a function of x and ip, solutions of equation 



p 1 

du = — dx H — <i-0 



exist provided 
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d (I 



dip \ qj dx \q J 
Eq. (|44] ) leads to the solvability condition 

2 dp 1 dg 2 ldg 2 = Q 
9^ 2 dip 2 dx 



(42) 



(43) 



(44) 



(45) 



Substituting g 2 and p from Eqs. fl3"T| ) and ( |4"2"D in Eq. ( [4"5] ) yields a rela- 
tion of the form Z]J=o a j{'4 , ) x ^ = 0; hence, since £ and are independent 
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variables, the equations a 3 - = (J = 0, . . . , 7) should be satisfied for all j. 
Consequently, the equilibrium problem is reduced to a set of eight ordinary 
non- linear differential equations for if) containing the surface functions P(if>), 
P s (if)), F(if)), &(if)), X(if)), p(if>) and r(if>). It is therefore over-determined. 
Furthermore, to completely solve the equilibrium problem, one should have 
to obtain the function u(x, if>) which by Eq. (fl3| ) satisfies the relation 



du 
dx 



P(x,if)) 
q(x,if;) 



(46) 



where p(x, if)) is given by Eq. ( 
q(x, if)) = xh 2 



and 



2(< + ^3 + ±m) 



(47) 



by Eq. (j37D. Therefore, the existence of T = T(if>) helically symmetric 
equilibria with incompressible flows, remains an open question. It may be 
noted that from the present and our previous studies |], |7j it turns out 
that, as concerns the mathematical structure and the physical properties of 
T = T(if>) equilibria with incompressible flows, there are similarities with 
static omnigenous equilibria. The non-existence of three dimensional static 
omnigenous equilibria has been proved in an elegant manner by Palumbo and 
Balzano f2(| ^TJ. In the presence of flow, however, there are difficulties in 
applying their method because of the departure of the isobaric and current 
density surfaces from the magnetic surfaces. 



IV. Conclusions 



The equilibrium of helically symmetric plasmas with incompressible flows 
has been investigated within the framework of ideal magnetohydrodynamic 
theory. For the system under consideration the equilibrium equations reduce 
to a second-order elliptic partial differential equation for the helical mag- 
netic flux function if) [Eq. (|25|)1 containing the density p(if>), the electrostatic 
potential $('?/'), the static equilibrium pressure -P S (V'), the function F(if>) as- 
sociated with the helical flow and the function X(if>) related to the magnetic 
field component along the symmetry direction in conjunction with a relation 
for the pressure [Eq. (0)]. This equation has been transformed to one, 
possessing differential part identical to the corresponding static equilibrium 
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equation [Eq. (B7T)] , which permits the construction of several classes of non- 
constant Mach-number analytic solutions. We have obtained exact solutions 
with flows in the direction of symmetry, parallel to the magnetic field and in 
arbitrary direction with differentially varying electric fields perpendicular to 
the magnetic surfaces. 

For the physically appealing case of equilibria with isothermal magnetic 
surfaces it has been shown that the equilibrium equations become over- 
determined, viz. the equilibrium problem reduces to a set of eight ordinary 
differential equations for if) with seven surface quantities. This result is dif- 
ferent from the corresponding ones of cylindrical and axisymmetric equilibria 
for which the condition T(ip) imposes only restrictions to the shapes of the 
magnetic surfaces. Thus, the existence of helically symmetric equilibria with 
isothermal magnetic surfaces as well as of generic three dimensional equilibria 
of this kind remain open questions. 
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